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ABSTRACT 

An understanding of icing physics is required for the de- 
velopment of both scaling methods and ice-accretion 
prediction codes. This paper gives an overview of our 
present understanding of the important physical proc- 
esses and the associated similarity parameters that de- 
termine the shape of Appendix C ice accretions. 

For many years it has been recognized that ice accretion 
processes depend on flow effects over the model, on 
droplet trajectories, on the rate of water collection and 
time of exposure, and, for glaze ice, on a heat balance. 
For scaling applications, equations describing these 
events have been based on analyses at the stagnation 
line of the model and have resulted in the identification 
of several non-dimensional similarity parameters. The 
parameters include the modified inertia parameter of the 
water drop, the accumulation parameter and the freezing 
fraction. Other parameters dealing with the leading- 
edge heat balance have also been used for conven- 
ience. By equating scale expressions for these parame- 
ters to the values to be simulated a set of equations is 
produced which can be solved for the scale test condi- 
tions. Studies in the past few years have shown that at 
least one parameter in addition to those mentioned 
above is needed to describe surface-water effects, and 
some of the traditional parameters may not be as signifi- 
cant as once thought. Insight into the importance of 
each parameter, and the physical processes it repre- 
sents, can be made by viewing whether ice shapes 
change, and the extent of the change, when each pa- 
rameter is varied. Experimental evidence is presented 
to establish the importance of each of the traditionally- 
used parameters and to identify the possible form of a 
new similarity parameter to be used for scaling. 

INTRODUCTION 

Many components cannot be tested full size in an icing 
wind tunnel because of test-section blockage limitations. 
Furthermore, facilities that simulate natural icing can 
provide only limited ranges of air speed, cloud drop size, 
and liquid-water content. To accommodate sub-scale 
models or to extend the useful range of testing, it is often 
necessary to scale the model size or test parameters. 
The objectives of such scaling are to insure that the 
amount of ice relative to the model size, the shape of the 
ice accreted, dimensionless ice accretion limits, and the 
aerodynamic penalties due to the ice are the same as 
would have been obtained with the desired, or refer- 
ence, model size or test conditions. It has always been 


assumed that for adequate scaling, characteristic fea- 
tures, such as glaze horns, need to be simulated in size, 
location, angle and shape, but how closely these charac- 
teristics need to match between scale and reference 
accretions has never been shown. For testing of ice- 
protection systems, the main objectives may be to match 
scale and reference non-dimensional accretion limits or 
water loading. In any case, it is important that scaling 
methods account for the appropriate physics. 

Accurate scaling requires that similitude of geometry, 
flow field, drop trajectory, drop catch and heat balance 
be satisfied. However, these relationships may some- 
times lead to conflicting values for one or more test pa- 
rameters, making it necessary to apply simplifications to 
permit their use in practical scaling situations. A number 
of scaling methods have been derived from the basic set 
of similitude equations with emphasis on different simi- 
larity parameters and with varying assumptions applied. 
These methods provide guidance for the researcher in 
establishing proper scaling conditions. 

For the testing of airfoils, when the desired model size 
cannot be accommodated in a test facility, an alternative 
approach to attempting to apply simplifications to the 
similitude equations is to use a hybrid scaling 
method . 123 In this approach, a full-size leading-edge 
portion of a wing section is attached to a flapped section 
of reduced size. The flap is adjusted such that the flow 
field around the leading-edge region of the scale model 
is the same as that of the full-size airfoil. In the hybrid 
method, scale test temperature, airspeed, liquid-water 
content, median volume drop diameter and icing time 
are the desired full-size values. However, multiple mod- 
els need to be fabricated for different flow conditions and 
angles of attack. The hybrid method will not be dis- 
cussed further here. 

In this paper, the most important similarity parameters 
will be reviewed, and results shown from icing tests per- 
formed to evaluate similarity parameters. The analyses 
and test results reported are for unprotected, unswept 
geometries for Appendix C conditions. Studies of the 
physics required to describe phenomena related to su- 
per-cooled large drop (SLD) icing is just beginning and 
will not be discussed. Nearly all of the experimental re- 
sults presented were obtained in sea-level tunnels. 

SIMILITUDE ANALYSIS 

A scaling method is a procedure to determine the scaled 
test conditions necessary to produce an ice shape that 
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simulates what the full-size model would accrete if ex- 
posed to the desired cloud conditions. To achieve this 
similarity in ice accretion, it is necessary for the scale test 
to simulate the geometry, the flowfield, the drop trajecto- 
ries, the total water catch, the surface heat balance and 
the surface phenomena of the desired icing encounter. 
Because of the potential complexity, the analysis in this 
presentation is limited to the stagnation line. The assump- 
tion is made that if the scale accretion here is similar to 
the reference case, so will that which forms everywhere 
on the model. 


This problem has usually been solved by scaling analy- 
ses by disregarding both of these parameters. The ra- 
tionale for this approach is that, for most icing condi- 
tions, the Mach number is relatively low and compressi- 
bility effects can be neglected. The Reynolds number 
has been ignored by arguing that ice accretions occur 
mostly at the leading edge of surfaces where the bound- 
ary layer is initially thin, and viscous effects are small. 
Farther downstream, as ice accretes it often causes the 
boundary layer to transition to turbulent flow, whose 
characteristics are nearly independent of Reynolds 
number. 


GEOMETRIC SIMILARITY 


The alternative to the hybrid scaling mentioned above is 
to make the reference and scaled models geometrically 
similar over the entire model. In this paper, the discus- 
sion will assume that scale and full-size, or reference, 
models have the same non-dimensional coordinates. As 
ice grows on the model, the shape must continue to be 
similar for both scale and reference to permit flowfield 
similarity. This basic requisite will be assumed for sub- 
sequent derivations of the other similitude requirements. 

FLOWFIELD SIMILARITY 


Similarity of the flowfield would suggest that the Rey- 
nolds and Mach numbers for the scaled test need to be 
matched to their respective reference values. The usual 
definitions for these parameters will be used: 


Re = 


Vdp a 


A a 


( 1 ) 


M = 


V 

RJ sl 


( 2 ) 


In equation (1) the length scale d represents either the 
diameter if the model is a cylinder or twice the leading- 
edge radius if the model is an airfoil. It was common 
practice in previous analyses to use the chord of the 
airfoil as the characteristic length; however, the main ice 
accretion occurs in the region near the leading edge. 
Therefore, in studying scaling similarity it seems reason- 
able to define similarity parameters in terms of a length 
representative of the leading-edge region. The conven- 
tion of using d to represent either cylinder diameter or 
twice the airfoil leading-edge radius will be used in the 
definition of all similarity parameters in this paper. 


Icing occurs in the atmosphere only within the limited 
range of absolute temperature, T st , of 420 to 492°R (233 
to 273 K). Thus, because yax\d R a are constant proper- 
ties of air, equating the Mach number for the scaled and 
reference cases requires that the velocities for the two 
be very nearly the same independent of scale model 
size. However, matching the Reynolds numbers re- 
quires the scale-to-reference velocity ratio to be ap- 
proximately proportional to the reference-to-scale model 
size ratio. Clearly, then, Mach and Reynolds numbers 
cannot be simultaneously matched when the scale 
model is other than full size. 


In addition to Mach and Reynolds number considera- 
tions, simulation of the flowfield requires that the veloc- 
ity, pressure and temperature distributions over the 
scale model must simulate those of the reference case. 
In subsequent derivations of equations, similarity of 
these distributions will be assumed. 

DROP TRAJECTORY SIMILARITY 

The mass of water reaching each part of the surface of 
each model must be similar. The drop trajectory deter- 
mines if an individual drop will impinge on the surface of 
the model and, if so, where on the body. Thus, similarity 
of drop trajectory between scale and reference cases 
should ideally be satisfied. An analysis of the drop tra- 
jectories will be made in this section, and expressions 
will be derived which relate the scaled and reference 
values of some of the operating and test parameters of 
importance. The quantity of water available to reach the 
surface will be discussed in the section Water Catch 
Similarity. 

Several simplifications are required to derive workable 
equations. It is assumed that the drop-size distribution 
in the cloud is the same for both scaled and reference 
cases and that the trajectories can be adequately repre- 
sented by that of one representative drop diameter. 
Bragg, et al 4 showed that for the typical distributions of 
drop size and relative velocity of icing experiments the 
drop motion and catch on a surface can be accurately 
represented by using the drop median volume diameter 
(MVD). Cloud drops involved in Appendix C icing are in 
the range of 10- to 50-|im diameter; therefore, the effect 
of gravity is small and can be neglected. The drop mo- 
mentum equation describing the motion of a single drop 
relative to the airstream has been given by both Bragg 5 
and Ruff. 6 


The modified inertia parameter, K 0 , was developed by 
Langmuir and Blodgett 7 by combining two dimensionless 
terms of the drop momentum equation: 


K n =^ + 



(3) 


where K is the drop inertia parameter: 


K = 


V5 2 Pw 

18 dp a 


(4) 
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and A/Astokes is the range parameter 


A 


A 


Stokes 


a Re s 

A_ f 

Re, l 


24 

^D^ e re/ 


dRe„ 


( 5 ) 


Here, 2 Stotes is the range of the drop if it were released in 
still air at an initial velocity of V with the drag given by 
Stokes' law, and A is the drop range with the actual drag. 
Re re , is the Reynolds number based on the relative ve- 
locity between the drop and air, and Re s is the drop Rey- 
nolds number: 


Re s 


V5p a 

Pa 


( 6 ) 


The range parameter defined by equation (5) was tabu- 
lated by Langmuir and Blodgett as a function of Re s . A 
fit to the Langmuir and Blodgett tabulation is 

i ( 0.8388 + 0.001 483Re s V 1 

— - — = , (7) 

^ Stokes [+0.1 847 Vr^ 

To satisfy drop trajectory similarity it is only necessary 
that K 0 ,s = K 0 ,r. Equation (3) is only valid if K is greater 
than 1/8. Langmuir and Blodgett noted that for K < 1/8 
drop impingement would not occur. 


Olsen 8 at NASA, Charpin and Fasso 9 at ONERA and 
others used a simplified form of equation (3) along with 
an approximation to permit integration of equation (5). 
Bragg 8 also presented alternative approaches to that 
given here. These simplifications permitted solving for 
the scale drop size directly in terms of reference condi- 
tions and other scale parameters. When equations are 
programmed into computers to solve for scaling condi- 
tions, however, these simplifications are unnecessary, 
and scaling studies at NASA Glenn for the past 10 years 
have matched scale and reference K 0 to solve for scale 
drop size. 


The catch efficiency, /?, varies over the model surface 
with a value p 0 at the stagnation line. Langmuir and 
Blodgett 7 published tables of the stagnation-point catch 
efficiency as a function of the parameters, K and K 0 for 
cylinders. They showed that fairly accurate values of p 0 
could be calculated from the following equation for K < 
7.5: 


1 .40 (K 0 - 1 / 8)° 84 
1 + 1 .40(K 0 — 1 / 8)° 84 


( 8 ) 


To test the validity of this expression for a range of icing 
conditions, p 0 calculated from equation (8) was com- 
pared with values from the LEWICE ice accretion code 10 
at both 0 and 10°AOA. The results are given in figure 1, 
and it can be seen that there is close agreement be- 
tween LEWICE and equation (8). For the NACA 0012 


airfoil used for these computations ji 0 varies little from 0 
to 10° AOA. If the scale model is mounted at the identi- 
cal AOA as the reference, matching of K 0 will still pro- 
duce the correct drop trajectories for the scale test al- 
though there will be a small error in p 0 if equation (8) is 
used. The conditions considered in figure 1 provide a 
range of inertia parameter, K, of 0.72 to 202, suggesting 
that Langmuir and Blodgett’s upper limit of K = 7.5 for 
equation (8) is very conservative. 

Equation (8) shows that to match p 0 between scale and 
reference values, it is only necessary to match K 0 . The 
initial capture efficiency is assumed to be that at the 
stagnation line of a clean surface. It can be argued for 
scaling, however, that aft of the stagnation line, because 
scale and reference models are geometrically similar 
and the flowfields are also similar, the collection efficien- 
cies must vary in the same way for both models. This 
agreement of scale and reference collection efficiencies 
is shown in figure 2, which gives LEWICE predictions for 
two NACA 0012 airfoils of different sizes. The condi- 
tions for the smaller model were scaled from those of the 
larger such that the stagnation K 0 (and therefore p 0 ) for 
the two matched. The two curves are indistinguishable 
over the range for which accretion occurs; to help iden- 
tify the curve for the reference case the area under that 
plot has been shaded. Figure 2 confirms that it is only 
necessary to match the stagnation K 0 to properly scale 
drop trajectories that impinge everywhere on the model. 
It is assumed that as ice accretes, because the geome- 
try changes in the same way for both scale and refer- 
ence models, the time-varying collection efficiency will 
continue to match everywhere. 



Modified Inertia Parameter, K Q 

Figure 1 . Stagnation Collection Efficiency for NACA 0012 
Airfoils at 0 and 10° AOA. t sh -30 to 25°F; p sf , 14 psia; V, 
100 and 400 mph; S, 10 to 100 pm; LWC, 1 g/m 3 . 
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Derivation of the equations for trajectory similitude as- 
sumed that a single drop diameter, the median volume 
diameter, adequately represents the cloud and its trajec- 
tory. At any instant during icing the size of a drop im- 
pacting a particular location on the model surface cannot 
be predicted. Thus, in glaze ice conditions in particular, 
features of the ice shape tend to be subject to a certain 
randomness; therefore, icing encounters repeated with 
the same model and identical test conditions cannot be 
expected to produce identical ice shapes. Fortunately, 
however, noticeable differences on the macroscopic 
scale are usually minor. 


Effect of Bn on Ice Shape 


In 1998 Chen 11 performed several icing tests to compare 
shapes with drop MVD’s of 20, 36 and 55pm. The tests 
were made in the IRT using a 24-in-chord GLC 305 air- 
foil model with constant t sh V and accretion time. LWC 
was adjusted to maintain constant freezing fraction, n 0 , 
for each set of tests. The freezing fraction was defined 
by Messinger 12 as the ratio of the mass of water that 
freezes at a given location on the surface to the total 
mass of liquid water that reaches the surface at that lo- 
cation. The freezing fraction will be discussed in more 
detail in the section on energy balance similarity. In 
Chen’s study, because velocity and model size were 
fixed, p 0 varied with MVD. The quantity of ice collected 
is determined by the product p 0 A c , and this value was 
maintained within about 10%. The accumulation pa- 
rameter, A c , will be discussed in the next section. Re- 
sults for drop sizes of 55 and 20 pm are given in 



-0.15 -0.10 -0.05 0.00 0.05 0.10 0.15 


Distance From Stagnation Line, s/c 

Figure 2. LEWICE-Predicted Reference and Scale Col- 
lection Efficiencies. NACA 0012 Airfoils. Reference 
Conditions: c R , 21 in (53.3 cm); V R , 150 mph (67 m/s); 
S R , 30 pm. Scale Conditions: c s , 10.5 in (26.7 cm); l/ s , 
262.5 mph (117 m/s); S s , 15.6 pm. 


figure 3. Although p 0 changed from about 92 to 74%, 
this reduction appeared to have no effect on the main 
ice shape. Undoubtedly, the icing limit would have 
changed with drop size (that is, with /? 0 ), but this feature 
was not measured. 

The implications of these results for scaling are signifi- 
cant. Within the FAA FAR-25 Appendix C envelope, the 
main ice shape appears to be virtually independent of 
drop size. Thus, if calculated conditions for a scale 
model test involve drop sizes outside the available range 
for a particular test facility, it may be possible to substi- 
tute a different 8, providing the key parameters PoA c and 
n 0 are maintained at the reference values. If p 0 is not 
matched, however, trajectories impinging aft of the main 
ice shape will not be similar for scale and reference. 
Therefore, when impingement limits are needed they 
must be found from tests whose conditions provide a 
match of p 0 . Further testing is needed to define better 
the effect of drop size on ice shape. 

WATER CATCH SIMILARITY 

The total amount of water reaching the surface depends 
on local collection efficiency, rate of water flow towards 
the model and accretion time. Assuming that all the im- 
pinging water freezes on impact (i.e., n = 1), that p is 
independent of the changing ice shape with time and 
that pi does not significantly vary as ice is accreted, the 
ice thickness, A would be given by: 

Pi 

At the stagnation line the ice thickness for glaze ice (n 0 < 
1) normalized by model size is given by 

y = AA" 0 (10) 
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S= 55 jim, Run 624.40 
5- 20 jim, Run 624.43 


Run 

in 

L St j 

°F 

v , 

mph 

47 , 

pm 

i— v v v_y, 

g/m 3 

h 

min 

% 

Ac 

n 0 

b 

°F 

17 , 

°F 

1 1C7, 
10 4 

10 6 

(a) 624.40 

24 

14 

200 

55 

1.16 

6.1 

91.6 

2.53 

0.28 

0.84 

16.1 

20.3 

10.9 

2.02 

624.43 

24 

14 

198 

20 

1.31 

6.1 

73.9 

2.83 

0.30 

0.76 

15.9 

20.2 

10.8 

1.98 

(b) 924.40 

24 

2 

200 

55 

1.16 

6.1 

91.7 

2.53 

0.51 

0.84 

28.6 

36.8 

11.4 

2.02 

924.43 

24 

2 

200 

20 

1.30 

6.1 

74.0 

2.83 

0.54 

0.76 

28.5 

36.7 

11.4 

2.00 


Figure 3. Effect of Drop MVD on Ice Shape. GLC 305 Airfoil. Ice-Shape Data from Tests by 
Chen. 11 


where the accumulation parameter, A c , is defined as 


LWCVt 


not always unity, and 
energy-balance simi- 
larity must be satis- 
fied for glaze ice ac- 
cretions. 


The rate at which 

water freezes on a 
— 1 1 1 1 1 1 1 1 1 1 surface depends on 

D.00 0.04 0.08 0.12 0.16 the local heat bal- 

ance. Tribus, et al 13 
studied heat transfer 

S= 55 urn, Run 924.40 for a heated c y |inder 

„ in icing conditions. 

- *= 20 Hm, Run 924.43 Their a " a|ysis i9 us0 . 

q g 2 ful for thermal deic- 

ing or anti-icing sys- 
<j>, 6, Re, We L , terns. Messinger 12 

J o D p iq 4 1q 6 later performed the 

8 0.84 16.1 20.3 10.9 2.02 classical energy- 

;0 0.76 15.9 20.2 10.8 1.98 , balanCe a " alySiS ,0r 

freezing at an un- 

1 0.84 28.6 36.8 11.4 2.02 heated surface. 

4 0.76 28.5 36.7 11.4 2.00 Messinger's work 

e-Shape Data from Tests by waj \ 6 de ^' 

K y oped by Ruff. The 

similitude relations to 

be presented here 

are for an unheated surface and are based primarily on 

the Messinger and Ruff analyses. The energy terms for 

this case, expressed per unit area, are: 


One goal of properly scaling the ice accretion is that the 
local normalized thickness of ice should be the same for 
the scaled and reference cases. If p 0 ,sA c ,sno,s = 
Po,rA c ,fPo.r the thickness of the scale ice at stagnation, 
relative to the model size, will match that of the refer- 
ence. We have already seen that to insure trajectory 
similarity it’s necessary to match scale and reference 
values of K 0 , and from equation (8) if the K 0 match, so 
do the p 0 - In the next section, it will be shown that n 0 
must also match. Therefore, water catch similarity (ice- 
thickness similarity) is typically satisfied by requiring that 
A cS = A c .fl. With the assumption that p, is the same for 
scaled and reference cases, the ratio of scale-to- 
reference accretion time from equation (11) is: 


fd s ) 

fM 

( lwc r \ 

7 

[Vs) 

{lwc s ) 


For airfoils, of course, the ratio ( ds/d R ) = ( Cs/Cr ). If it is 
not possible to find scaled conditions that permitted a 
match of K 0 , then scale and reference j3oA c must match. 

ENERGY BALANCE SIMILARITY 

If all the water impinging on a surface were to freeze on 
impact, K 0 ,s = K 0 , R and A CiS = A C , R would be sufficient to 
satisfy scaling. This situation applies to rime ice forma- 
tion. However, the fraction of water that freezes, n, is 


(1) Energy lost from the surface by convection through 
the boundary layer: q c = h c (t s - t b j). 

(2) Energy lost from the surface due to evaporation of 
water: q e = m e A v . 

(3) Energy lost from the surface to raise the tempera- 
ture of the impinging liquid to the freezing point: 

(4) Energy gained by the surface due to release of la- 
tent heat of fusion from the freezing water: q f = mn 0 A f . 

(5) Energy gained by the surface from the kinetic en- 
ergy of the water drops striking the surface: 

When the surface temperature reaches steady state, the 
net heat transfer is zero. Thus, the energy-balance 
equation is 

q c +q e +q w =q f +q k (13) 

As with the drop trajectory analysis, the energy balance 
will be evaluated only at the stagnation line for scaling. 
It is assumed that if scale and reference conditions are 
matched there, they will also match everywhere. In ad- 
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dition, the energy balance on the clean airfoil (before ice 
begins to accrete) is used. If the clean heat balance 
matches for two models and similarity is maintained in 
all the physical phenomena that affect the shape of the 
accreting ice, the balance, although changing with time, 
will be the same for the two models as ice develops. 

Equation (13) is adequate for a scaling analysis. Other 
terms needed for a complete analysis include the sensi- 
ble heat gained as ice cools from the freezing tempera- 
ture to the surface temperature, sublimation, radiation, 
surface water in-flow, surface water run back and con- 
duction into the model surface. For glaze ice there is no 
sensible heat from cooling the ice, because the surface 
temperature, f s , is equal to the freezing temperature, t f . 
For this paper, the energy balance will be written for 
glaze ice only, so this contribution is neglected. The 
sublimation is small compared with the evaporation of 
liquid, and for unheated surfaces radiation is also insig- 
nificant. The analysis is limited to the stagnation line of 
the model where water arrives only by impingement; 
therefore, heat carried into the control volume by water 
flowing from a neighboring location can be omitted. 
Heat carried away by water flowing out of the control 
volume to a neighboring location has also been ne- 
glected on the basis that temperature differences be- 
tween a stagnation control volume and a neighboring 
one will be negligible. For thermal ice protection this 
term may need to be considered, but for this paper ice- 
protection systems will not be discussed. Finally, con- 
duction effects can be neglected for icing times greater 
than a few seconds because the ice quickly produces an 
insulating layer over the model surface. 

The film coefficient in term (1), h c , varies with position on 
the model. Kreith 14 gives the following relationship for 
Nu for the stagnation line of a cylinder: 

Nu = / \A4Pr 0A Re 05 (14) 

where 

.. h d 

Nu=-^~ (15) 

k 


validity was confirmed by Poinsatte 15 and Van Fossen. 16 
At the stagnation line of a clean, smooth NACA 0012 
airfoil at 0° AOA in the IRT, Poinsatte found 

Nu = I.IOORe 0 472 (18) 

The coefficients in Poinsatte’s study had an uncertainty 
of about ±5%. Poinsatte’s original results used Nu and 
Re based on chord, but in equation (18) they have been 
converted to a c/-basis for consistency with the practice 
in this paper. For the NACA 0012 airfoil, d = .0316c. 
The Poinsatte 0012 data correlations of equation (18) 
are within 10% of the Kreith cylinder coefficients and 
exponents of equation (17). Van Fossen’s clean, 
smooth cylinder tests for IRT turbulence levels produced 
a coefficient about 25% higher than Kreith’s and a Re 
exponent within 5%. For scale testing in an icing tunnel 
with turbulence levels similar to those of the IRT, it is 
recommended that Kreith’s expression for Nu be used 
because of its consistency with these experimental re- 
sults for both cylinder and airfoil models. 

At the stagnation point, t u = t st + \fl2c P:a . Thus, heat 
balance term (1 ) is q c = h c (t s - t s , - \fl2c p J). 

In heat balance term (2), the mass of water which 
evaporates, rh e , is 

m e = h J Pww ~ Pw ^ (19) 

Pst 

The driving force for evaporation in equation (19) is the 
difference between the partial pressure of vapor at the 
surface, p ww , and that in the atmosphere, p w . The form 
given here neglects compressibility effects. The vapor 
pressure over water at the surface is assumed the satu- 
ration pressure at the surface temperature. A curve fit 
for the saturation pressure of vapor over water was 
given by Pruppacher and Klett 17 for the range of -50°C to 
50°C. Their expression is a convenient way to evaluate 
p ww and p w . Ruff 6 derived a compressible form of the 
evaporation driving force. Equation (19) in compressible 
form is 


The air properties for Nu, Pr, and Re in equation (14) are 
evaluated at the film temperature, which is taken as the 
average of the free-stream and surface temperatures, 

t, m =l(t s + t sl ) (16) 

For glaze ice the surface temperature is 0°C, and over a 
free-stream temperature range of -40 to 0°C, the range 
of interest to icing, the Pr evaluated at film temperatures 
varies only from .708 to .705; thus, 1.14Pr 4 has a nearly 
constant value of .992. Then 

Nu = 0.992Re° 5 (17) 

Ruff 6 used equation (14) to find the convective heat 
transfer coefficient for both cylinders and airfoils and this 
expression is also used in NASA scaling research. Its 



Compressibility effects are typically negligible for icing 
conditions. In equations (19) and (20) h G is the gas- 
phase convective mass transfer coefficient. It is analo- 
gous to the convective heat-transfer coefficient. The 
theoretical relationship between h G and h c can be ex- 
pressed in terms of the Schmidt and Prandtl numbers as 
(see, for example, ref. 14): 



( 21 ) 
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In equation (21) Pr and Sc are those for air, with proper- 
ties evaluated at the film temperature, from equation 
(16): 


p r C P.aP a 

k a 

(22) 

Sc = ^ a 
P,D„ 

(23) 


With these definitions, the energy balance of equation 
(1 3) written for the stagnation line is 


L -f„, - 


V 2 

2c 


1~tot Pst 


P.a J 


V 


1 P,o 
0.622 T„ 


A, 


(24) 


St J 


, . .v 2 

+ i™ p .ws ( f f - ) = mnA +rn— 


Tribus 13 introduced the dimensionless parameter b 
known as the relative heat factor: 


mc pws LWCV f3 0 c p , 
ho K 


(25) 


It relates the total heat capacity of the impinging water to 
the ability of the surface to convect heat. Note that b is 
evaluated at the stagnation line because J3 0 is the local 
stagnation catch efficiency, found from equation (8). 


Two other parameters often used for convenience are <j) 
and 6 which have dimensions of temperature and relate 
to the drop energy transfer and air energy transfer, re- 
spectively: 


and 


(P = t f -t s , 


V 2 


(26) 


0 = 


t - r - 


v 2 

2c 


T tot Pst 


P.a J 


1 


A, (27) 


0.622 T 


St J 


This form of # includes compressibility effects. 


When the glaze-ice energy balance, equation (24), is 
rewritten in terms of the parameters b, (/and 0, the terms 
can be arranged to give an expression for the freezing 
fraction: 


n 


0 


f C P.«s ] 


A, 



(28) 


A freezing fraction of 0 implies that no freezing occurs, 
while a freezing fraction of 1 indicate that water freezes 
on impact. The latter situation produces rime ice. Ex- 
perimental evidence to validate equation (28) will be pre- 
sented below. 


The Ruff scaling method finds scale LWC by setting n 0 ,s 
= n 0 ,R. By using n 0 as the heat-balance parameter to 
match for scale studies, the implicit assumption is that if 
n 0 ,s = no, a then n will also match everywhere. Further- 
more, strictly speaking, equation (28) applies only to a 
clean airfoil. As ice accretes, the stagnation freezing 
fraction can be expected to change because the geome- 
try is continually changing. Thus, for scaling it is neces- 
sary to assume additionally that if n 0 varies with time, it 
does so in such a way that if n 0 ,s = n o.R at the initiation of 
spray, freezing fractions will continue to match during 
the course of ice accretion. This is not an unreasonable 
assumption because the goal of scaling is to maintain 
similar ice shapes for scale and reference. If this goal is 
met the time variation of n 0 ,s and n 0 ,R should be the 
same. 

Kind and coauthors' 18 ' 19 formulated the energy equation 
to give a freezing-fraction expression valid for any posi- 
tion on the surface, rather than just the stagnation line. 
Two parameters were used with the first incorporating 
terms independent of location and the second varying 
with distance from the stagnation line. They gave 
LEWICE predictions showing that when the stagnation 
n 0 matched, reference and scale local freezing fractions 
also matched from the stagnation line back to about s/c 
= ±0.02. This is the region of greatest interest to icing 
because most of the main ice shape is accreted there. 
However, it’s important to keep in mind that LEWICE 
predictions of freezing fraction have not been validated 
experimentally the way collection efficiencies have, and 
therefore there is less confidence in the value of n from 
LEWICE. 

Kind and associates also concluded that scale ice 
shapes simulated reference shapes best when the scale 
and reference static temperatures matched. However, 
all of their observations were made with scale velocities 
of 100 m/s (224 mph) or less. For such air speeds the 
velocity term is never more than approximately 2°F, and 
equation (26) shows that if scale and reference </> match 
the temperatures will also nearly match. It would appear 
that the matching of (/> is more general than simply 
matching t st . Matching of </> to find scale static tempera- 
ture is a requirement of the Ruff 6 scaling method. 

For rime ice, because water freezes on impact, the en- 
ergy balance does not affect the ice shape. Conse- 
quently, no energy balance needs to be considered 
when scaling rime. 

Effect of Temperature and LWC 

The strong effect of the freezing fraction on ice shape 
has been demonstrated by numerous studies with both 
cylinders and NACA 0012 airfoils 20 ' 21 . A few examples 
of the effect of changing either temperature or LWC, 
both of which have a significant influence on n 0 , will be 
given. 

Ice shape changes with incremental changes in tem- 
perature are shown in figure 4. These unpublished tests 
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were performed in the NASA Glenn IRT using 21 -in- 
(53.3-cm-) chord NACA 0012 airfoil models. In each 
portion of the figure, two ice shapes are compared for 
which the temperature was changed with all other test 
parameters held constant. For all tests K 0 was constant 
to give a fixed value of j8 0 of 0.85, A c was 1 .88 and b was 
0.58. For each 5°F decrease in temperature, the pa- 
rameter </> increased approximately 5°F and 6 increased 
about 7°F. The freezing fraction also increased as indi- 
cated in the figure with decreasing temperature. The ice 
shape for the warmer-temperature of each pair is shown 
shaded while the colder of the two is represented with a 


With each incremental decrease in temperature (and 
consequent increase in n 0 ) the horns moved forward and 
the leading-edge ice thickness increased (consistent 
with eq. (10)). The exception to this general trend oc- 
curred when n 0 was increased from 0.52 to 0.61 (fig. 4 
(c)). Although the ice thickness increased, there was no 
significant change in horn angle for this range of freezing 
fraction. 

With constant test parameters c, V, 8, LWC, and rand 
constant similarity parameters K 0 (or p 0 ), A c and b for the 
range of temperatures tested in figure 4, it must be rea- 
soned that only the varying parameters n 0 , </>or #can be 

responsible 



for the 
changes 
observed in 
ice shape. 
Further evi- 
dence will 
next be 
shown to 
isolate 
which of 
these three 
has the 


-0.04 0.00 0.04 0.08 -0.04 0.00 0.04 0.08 -0.04 0.00 0.04 0.08 greatest 


x/c 


x/c 


impact on 
shape. 


t s( , 19°F (-7°C); n 0 , 0.28 
f s(> 14°F(-10°C); n 0 , 0.40 

(a) Temperatures of 19 and 14°F 
(-7 and -10°C). 


f s( , 14°F (-10°C); n 0 , 0.40 t st , 9°F <-13°C); n 0 , 0.52 

t sl , 9°F (-13°C); n 0 , 0.52 f sf> 5°F <-15°C); n 0 , 0.61 

(b) Temperatures of 14 and (c) Temperatures of 9 and 5°F 
9°F (-10 and -13°C). (-13 and -15°C). 


Figure 5 

presents 
results of a 
series of 


solid line. 

0.08 
0.04 
0.00 
-0.04 
-0.08 

-0.04 0.00 0.04 0.08 0.12 



tests re- 
ported in reference 20. For these tests, the parameters 
c, V, and 8 were maintained constant; consequently, K 0 



-0.04 0.00 0.04 0.08 0.12 


(or po) was constant, while LWC 
varied. The icing time was ad- 
justed to keep A c the same for 
both tests in each portion of the 
figure. Again, the model was a 
21 -in- (53.3-cm-) chord NACA 
0012 tested in the NASA Glenn 
IRT. The shaded ice shape in 
each portion of the figure was 
recorded for an LWC of 1 .4 g/m 3 , 
while the shape indicated by a 
solid line was obtained with an 
LWC of 0.8 g/m 3 . 


x/c x/c 

t sl , 5°F (-15°C); n 0 , 0.61 t st , -2°F (-19°C); n 0 , 0.77 

t s! , -2°F (-19°C); n 0 , 0.77 t st , -16°F (-27°C); n 0 , 1.00 

(d) Temperatures of 5 and -2°F (e) Temperatures of -2 and -16°F 

(-1 5 and -1 9°C) (-1 9 and -27°C). 

Figure 4. Effect of Static Temperature on Ice Shape. Unpublished NASA data by 
Anderson. NACA 0012 Airfoil, c, 21 in (53.3 cm); V, 150 mph (67 m/s); 8, 38 pm; 
LWC, 1.0 g/m 3 ; r, 7.3 min; p 0 , 0.85; A c , 1.88; b, 0.58. 


For the two ice shapes shown in 
figure 5 (a) t s , was the same; 
therefore, n 0 increased with de- 
creasing LWC. As the LWC was 
decreased from 1.4 to 0.8 g/m 3 
the horns on the main ice shape 
moved forward and the leading- 
edge thickness increased, similar 
to the effect of decreasing the 
static temperature shown in 
figure 4. For both shapes $ was 
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20°F and <9 was 28°F. The relative heat factor was 0.77 
for the higher LWC test and 0.44 for the lower LWC. 
Therefore, keeping <j> and 0 constant failed to keep the 
ice shapes the same when freezing fraction and b var- 
ied. 

Figure 5 (b) gives results of tests in which the tempera- 
ture was varied to maintain a constant freezing fraction 
as LWC changed. These two ice shapes agreed. The 
value of (j> changed from 27 to 18°F and 6 decreased 
from 36 to 25°F with the increase in temperature from 4 
to 13°F. b decreased from 0.77 to 0.44 with the reduc- 
tion in LWC from 1.4 to 0.8 g/m 3 . The ability to match 
shapes by matching n 0 , although the temperature and 
liquid water content differed, shows that t st and LWC do 
not have an effect on ice shape independent of n 0 . Fur- 
thermore, as long as freezing fraction matches, it is not 
necessary to match the energy parameters </>, dor b. 

The results of the three test series shown in figures 4 
and 5 demonstrate the very strong effect of freezing 
fraction on ice shape. The energy parameters b , <^and 9 
are related to n 0 through equation (28) but appear to 
have little or no independent effect on ice shape. These 
results demonstrate that either LWC or t sl can be scaled 
from one value to another simply by adjusting t st or LWC, 
respectively, to keep the freezing fraction constant. This 
approach to scaling test conditions is known as the Ol- 
sen method. 8 

Validation of the Freezing Fraction Expression 

Using equation (10) it is possible to calculate an experi- 
mental freezing fraction from the leading-edge thickness 
once the B 0 and A c are determined from test conditions. 
This was done by Anderson and Tsao and the results 
compared with the analytical freezing fraction formula- 
tion presented in equation (28). To distinguish between 
the experimentally-derived and analytical freezing frac- 
tions, n e and n a will be used, respectively. 


The Anderson and Tsao study was based on icing tests 
in the IRT using NACA 0012 models with chords of 10.5 
- 31 .5 in (26.7 - 80 cm), velocities of 1 06 - 257 mph (48 
-115 m/s) and drop MVD’s of 21.5 - 46 pm. The ex- 
perimentally-derived freezing fractions are plotted 
against the analytical values in figure 6. Open symbols 
were measured from ice tracings taken at mid-span (CL) 
and the solid from the tracings 1 in (2.5 cm) above mid- 
span (CL+2.5 cm). The solid line represents perfect 
agreement of the two freezing fractions and the shaded 
band indicates the limits for ±10% agreement. A linear 
fit to the data is shown as a dashed line. 

The linear fit to the data fell within the ±10% agreement 
band, although on average n e values tended to be about 
5% lower than the n a . The value of n e found from equa- 
tion (10) is inversely proportional to the value of LWC 
through the accumulation parameter (eq. (11)). Be- 
cause the uncertainties for the LWC and A c were both 
estimated to be about ±10%, this agreement is quite 
good. 

Bilanin 23 performed a similar study using a limited set of 
ice-shape data from Ruff. Bilanin found agreement of n e 
and n a at rime, but the two deviated significantly for the 
lowest freezing fractions. He concluded that the Mess- 
inger analysis for freezing fraction had serious deficien- 
cies, probably because surface-water effects, particu- 
larly splashing, were not included in the heat balance. 
However, Bilanin did not consider the possibility that the 
facility LWC calibration or the uncertainty in the leading- 
edge thickness for low freezing fractions might have af- 
fected his results. 

Unlike Bilanin’s results, the experimentally-determined 
freezing fractions for the Anderson and Tsao study in 
figure 6 showed no systematic deviation from the ana- 
lytical values at low freezing fractions. This consistently 
good agreement between n e and n a over a range of con- 
ditions including both fully rime and fairly warm glaze 
argued against the contention that the Messinger freez- 



LWC = 1.4 g/m 3 , n 0 , 0.40 LWC = 1.4 g/m 3 , f s( , 4°F 

LWC = 0.8 g/m 3 , n 0 , 0.59 LWC = 0.8 g/m 3 , t sP 13°F 


(a) Constant Static Temperature. (b) Constant Freezing Fraction. 

f sf =11°F(-12°C) n 0 = 0.52. 

Figure 5. Effect of LWCor\ Ice Shape. 20 NACA 0012 Airfoil; c, 21 in (53.3 cm); 
V, 150 mph (67 m/s); S, 30 pm; r, 5.2 - 7.3 min; p 0 , 0.80; A c , 1.90. 
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Freezing Fraction from Analysis, n a 


Figure 6. Experimental and Analytical Freezing 
Fractions of Anderson and Tsao. 22 

ing fraction formulation neglects any important surface- 
water effects. On the contrary, for Appendix C icing 
conditions, the experimental results appeared to validate 
the Messinger freezing fraction formulation expressed in 
equation (28). 

Besides the consistency between n e and n a over the 
range of freezing fractions evaluated, Anderson and 
Tsao showed that when conditions were such that equa- 
tion (28) gave n a = 1 , ice shapes as well as the white 
appearance of ice were typical of fully rime ice. The 
rime shape in figure 7 is shown with a shaded cross sec- 
tion. A second shape, produced at a temperature giving 
an analytical freezing fraction of 0.91, is represented 
with a solid line. This shape differed from the rime in 
that a narrow strip of ice at the leading edge changed 
from the characteristic white appearance of rime to a 
slightly transparent form. At the same time the smooth, 
convex shape of rime was replaced at the leading edge 
with a slight valley. The remainder of the ice aft of this 
narrow strip maintained the rime appearance. 

Thus, the use of the Messinger analysis leads to a cal- 
culated freezing fraction of 1 when ice shapes are in- 
deed fully rime and less than 1 when glaze features be- 
gin to show at the leading edge. These results again 
show the strong consistency between the Messinger 
freezing fraction of equation (28) and experimental ice 
shapes. 

Anderson and Tsao also looked at the effects of calcu- 
lating the freezing fraction with 2 minor modifications. 
The first was to use the model chord instead of twice the 
leading-edge radius, and the second substituted the av- 
erage Nu around a cylinder 14 


0.08 i 


0.04 - 
0.00 - 
-0.04 - 


-0.08 H — i — i — i — i — i — i — i — i — i — i — i — i 

-0.08 -0.04 0.00 0.04 0.08 0.12 0.16 

x/c 

r> a =1 .00, 10-12-02 Run 11 
n 3 = 0.91, 10-12-02 Run 8 

Figure 7. Effect of Freezing Fraction on Ice Shape at n 0 
near 1 . NACA 001 2 Airfoil; c, 21 in (53.3 cm); V, 1 1 5 
mph (51 m/s); 8, 40pm; LWC, 1 .14 g/m 3 ; t st , -13 and -7°F 
(-25. and -22°C); J3 0 , 0.84; A c , 2.43; b, 0.57; </>, 40 - 44°F 
(21 - 25°C); 8, 52 - 59°F (29 - 33°C). 



Nu = 0.0239Re° 8OS (29) 

for the stagnation value of equation (14). Both of these 
changes led to somewhat poorer agreement between 
the values of n e and n a , but more significantly gave ana- 
lytical freezing fractions of 0.82 to 0.89 for the icing tests 
that produced the fully rime ice shown in figure 7. 

These comparisons illustrate how it is possible to apply 
simple checks of experimental consistency to expose 
problems with particular analytical methods for calculat- 
ing freezing fraction. At the same time, because fully 
rime is easy to recognize, equation (10) provides a sim- 
ple way to validate facility LWC calibrations. For this 
purpose, tests should be made at temperatures cold 
enough to result in fully rime ice. With /3 0 determined 
from equation (8) and the stagnation ice thickness 
measured from the experimental accretion, equation 
(10) can be solved for the actual A c , and the true LWC 
found from equation (1 1 ). 

SIMILARITY OF SURFACE-WATER DYNAMICS 

Scaling methods developed from the 1950’s through the 
1970’s relied on matching scale and reference values of 
some or all of the parameters K 0 , (or fi 0 ) A c , and n 0 . 
Given a set of reference conditions and choosing scale 
model size, Ruff used these parameters plus ^ and 0 
from equations (26) and (27) to provide five equations to 
solve for scale temperature, pressure, drop MVD, cloud 
LWC and spray duration. All of these earlier methods, 
including Ruffs, left the arbitrary choice of scale velocity 
to the user. In this section, it will be shown that velocity 
has an effect on ice shape independent of its effect on 
the similarity parameters identified so far. This effect of 
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velocity appears to be associated with surface-water 
phenomena, as will be seen below. 

In 1988 Bilanin 23 presented a Buckingham-n analysis in 
which he concluded that surface-water phenomena had 
to be included in icing scaling methods. Olsen and 
Walker 24 and Hansman, et al, 26,27 studied surface ef- 
fects and surface water during ice accretion, presenting 
additional evidence that these were important phenom- 
ena to consider in ice accretion. From close-up photo- 
graphic studies, it was observed that for glaze ice accre- 
tion unfrozen water on the ice surface tended to coa- 
lesce to form beads. These beads sometimes were 
swept downstream and sometimes froze in place. Bi- 
lanin 23 ' 28 also argued that drop splashing on impact 
might affect the shape of the ice accreted. 

Hansman and Turnock 26 found that when a surfactant 
was added to the icing spray water, the ice appearance 
and shape changed significantly, with the glaze horns 
moving toward the leading edge. Similar observations 
were made by Bilanin and Anderson. Figure 8 includes 
published 28 and unpublished data from their studies and 
compares ice shapes with and without surfactant addi- 
tion to the spray for velocities of 105, 150 and 209 mph. 
The shaded ice shapes were obtained with the NASA 
Glenn IRT de-mineralized water supply with no surfac- 


tant added, and the shape represented by the solid line 
resulted from surfactant addition to the water. Other 
than the addition of surfactant the same test conditions 
were used for each pair of tests, as shown in the ac- 
companying table. Surfactant addition reduced the sur- 
face tension to roughly half that of water. Because f} 0 , 
A c and n 0 were the same for each pair of tests, the lead- 
ing-edge ice thickness was also nearly the same. How- 
ever the included angle between the horns decreased 
dramatically when surfactant was added. Horn angle 
also decreased when velocity increased although the 
freezing fraction was approximately the same for all 
tests. Thus, both surface tension and velocity have an 
effect on ice shape independent of freezing fraction, and 
scale velocity cannot be chosen arbitrarily. Clearly, 
then, a similarity parameter dependent on the ratio V 
e1 /a wa e2 must be included in scaling methodology, where 
the powers el and e2 still need to be determined. 

Kind, et al 29 ' 30 considered the capillary number, 
IJ W V lo wa , to describe the phenomena related to the sur- 
face water on glaze ice accretions. The capillary num- 
ber is the ratio of viscous forces to surface-tension 
forces. When scale tests use water with viscosity and 
surface tension that match the reference values, match- 
ing of the capillary number is equivalent to matching the 



x/c 


x/c 


x/c 


Water, 6-9-94 Run 4 


Water + Surfactant, 
8-30-93 Run 4 


(a) V = 1 05 mph (47 m/s) 


Water, 8-27-93 Run 9 

Water + Surfactant, 

8-30-93 Run 5 

(b) V= 150 mph (67 m/s) 


Water, 8-27-93 Run 8 
Water + Surfactant, 


8-30-93 Run 3 
(c) V = 209 mph (94 


Date/Run 

d, 

in 

tst, 

°F 

v ; 

mph 

s, 

pm 

LWC, 

g/m 3 

T, 

min 

£>wai 

dyne/ 

cm 

Po > 
% 

A c 

n 0 

b 

* 

°F 

e, 

°F 

Re, 

10 4 

We L , 

10 6 

(a) 6-9-94/4 

2 

18 

105 

40.0 

1.16 

16.0 

65 

66.0 

1.13 

0.28 

0.76 

13.7 

20.3 

18.5 

1.73 

8-30-93/4 

2 

18 

105 

40.0 

1.17 

16.0 

30 

66.0 

1.13 

0.29 

0.76 

13.9 

20.7 

18.5 

3.71 

(b) 8-27-93/9 

2 

11 

150 

34.1 

1.39 

10.2 

65 

65.4 

1.22 

0.33 

1.08 

20.3 

28.3 

26.7 

3.51 

8-30-93/5 

2 

11 

150 

34.1 

1.39 

10.2 

30 

65.5 

1.23 

0.32 

1.08 

20.3 

28.1 

26.8 

7.66 

(c) 8-27-93/8 

2 

10 

209 

30.0 

1.10 

9.0 

65 

65.6 

1.20 

0.31 

1.03 

20.2 

25.2 

36.6 

6.94 

8-30-93/3 

2 

10 

209 

30.0 

1.10 

9.2 

30 

65.6 

1.22 

0.31 

1.03 

19.6 

24.7 

36.4 

14.88 


Figure 8. Effect of Surfactant and Velocity on Ice Shape. Vertical Cylinders Tested in the NASA Glenn 
IRT. Published 28 and Unpublished Ice-Shape Data from Tests by Bilanin and Anderson. 
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scale and reference 
velocity. However, 
studies of the effect of 
scale velocity 31,32 
showed that for sub- 
scale models the 
scale velocity must be 
higher than the refer- 
ence. Therefore it is 
unlikely that the capil- 
lary number will prove 
to be an important 
similarity parameter 
by itself. 



x/c 

c- 36 in, Run 2-8-36 
c= 12 in, Run 2-8-12 


x/c 

c = 36 in, Run 0-14-8-3 
c = 12 in, Run 0-14-8-0 


(a) n = 0.4 (b) n = 0.5 
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(a) 

2-8-36 

36 

9.5 

265 

44 

0.58 

14.6 

87.2 

2.67 

0.40 

0.57 

19.5 

21.2 

21.4 

5.28 


2-8-12 

12 

8.5 

280 

26 

0.91 

3.0 

90.0 

2.68 

0.41 

0.55 

20.2 

21.2 

7.5 

1.96 

(b) 

0-14-8-3 

36 

4.9 

271 

38 

0.55 

13.5 

84.9 

2.40 

0.52 

0.54 

24.0 

26.7 

22.2 

5.54 


0-14-8-0 

12 

5.4 

270 

20 

0.96 

2.6 

85.9 

2.42 

0.50 

0.54 

23.5 

26.1 

7.4 

1.83 


Figure 9. Effect of Model Size on Ice Shape. GLC 305 Airfoil. Ice-Shape Data from Tests 
by Chen. 11 


Chen’s 11 results dis- 
cussed above (see 
fig. 3), and others like 
them for NACA 0012 
airfoils, show that 
drop size and collec- 
tion efficiency appear 
to have little effect on 
the main ice shape. 

Therefore, the pa- 
rameter being sought 
cannot be dependent 
on 8. Chen also 
evaluated the effect of model size, and figure 9 shows 
that with p 0 , A c , n 0 , b , ^ and 9 constant a reduction in 
chord from 36 to 12 in moved the glaze horns rearward. 
This is the same effect shown in figure 8 for decreasing 
velocity. 

The previous section showed that the temperature and 
LWC do not have effects on the ice shape independent 
of the freezing fraction. Thus, the supplementary pa- 
rameter does not need to include these test conditions, 
and its general form must be 

V e 'c e2 

P = constant — (30) 

® wa 

This form suggests a Weber number based on chord: 


Studies by Bartlett 33, 34 and Oleskiw, et. al. 35 found no 
measurable effect of pressure on ice shape. These ob- 
servations rule out the dependence on air density in 
equation (31), making water density a more likely 
choice. Furthermore, the length may not be chord itself 
but rather some physical characteristic related to the 
accreting ice that is proportional to chord; for example, 
the water-film thickness. Because this length is not yet 
identified, L will be used to represent it, as follows 

V 2 Lp 

We L = ^ (32) 

cr 


The results of figure 9 suggest that L <x c, or, since the 
leading-edge radius is proportional to chord, L oc d. 
Then, from equation (32) the scale velocity found from 
matching We LiS to We LiR is 



Some preliminary scaling results with We L matched 
were reported by Anderson and Tsao. 36 An example is 
shown in figure 10 for scaling from 36- to 21 -in chord for 
NACA 0012 airfoils. Note that the stagnation collection 
efficiencies for the two tests do not match, so the condi- 
tions do not strictly satisfy usual scaling requirements. 
However, it was shown in figure 3 that p 0 does not have 
a strong effect on the main ice shape. More important 
than p 0 matching is that the product PoA c was nearly the 
same for the tests shown. The shapes of figure 10 as 
well as others in reference 36 for the same scale ratio at 
freezing fractions of 0.3 and 0.4 showed good agree- 
ment between scale and reference. To date there have 
been no results published for greater scale ratios for 
which the We L was matched, and such data are needed 
for more confidence in this approach. 

Although the conditions for the tests of figure 10 left the 
(/> mismatched, the good agreement of the ice shapes 
suggests again that may not be as important a similar- 
ity parameter as once believed. Note also that the tem- 
peratures for the two tests were significantly different; 
this result gives further evidence that static temperature 
does not need to be matched for good scaling. 


NASA/TM— 2005-213851 


12 



0.08 n 



-0.08 i 1 1 1 1 1 1 1 1 \ \ 1 

-0.04 0.00 0.04 0.08 0.12 0.16 


x/c 

c, 36 in; 2-8-02 Run 7 
c, 21 in; 11-13-00 Run 5 

c, t sl , V, 8 LWC, t, Po, , 8 We L , 

in °F mph jam g/m 3 min % c 0 °F io 6 

36 -3 115 160 1.50 9.7 95.4 1.70 0.50 34 1.18 
21 9 150 38 0.99 7.3 84.9 1.88 0.52 22 1.16 

Figure 10 . Scaling Using We LS = We LR to Determine 
Scale Velocity. 36 NACA 0012 Airfoils. 

There are a variety of Weber numbers based on other 
lengths and different velocities which have been 

matched to find scale velocity. One such parameter, for 
example, 30,31,32 was based on drop size and water prop- 
erties, 

We = V 5Pw (34) 

®wa 

A form of equation (32) in which h fi i m , the water-film 

thickness, was substituted for L was evaluated by 

Anderson and Feo. 37 Nearly the same Weber number 
was proposed and evaluated by Kind and col- 
leagues, 18 ' 19,38 ' 39 but with p a used in place of p w . A We- 
ber number based on the velocity of the water film in- 
stead of the free stream could also be considered. 

In addition to these Weber numbers, Feo 40 suggested 
that the water film thickness might itself be a critical di- 
mension controlling surface water dynamics. From stud- 
ies of sprays impinging on a cylindrical probe at above- 
freezing temperatures he found that water-film thickness 
was proportional to the product LWC vz Re' V4 . 

A scale velocity obtained as the average of that found 
from matching We (eq. (34)) and that matching Re (eq. 
(1)) has also been considered. 31,32 For half-size scaling, 
the average velocity method gives a scale velocity that is 
only about 20% higher than the one found by matching 
We L . The use of other Weber numbers or of the water- 
film thickness nearly all give scale velocities within ±20% 
of that found from matching We L , even for scale sizes as 
small as 14 the reference. Choosing the best among the 
different parameters by comparing scale to reference 
shapes is not productive because ice shapes are not 


sensitive enough to velocity to be able to discriminate 
between these differences. Therefore, until additional 
research is done, it is not possible to positively identify 
the similarity parameter, or parameters, that best repre- 
sents surface physics, and finding scale velocity by us- 
ing equation (33) appears to be a satisfactory approach 
at the present time. 


OTHER PHENOMENA 

The physical models for icing on which the similitude 
parameters have been based have not been fully veri- 
fied and probably do not include all phenomena that may 
have some effect on the ice accretion process. Close- 
up studies of water impact and freezing on surfaces 
were made by Olsen and Walker. 24 They believed that 
shedding of water from the surface, particularly from the 
tips of the horns, is a significant part of the ice-accretion 
process. Whether this water is subsequently re- 
entrained and deposited aft of the main ice shape is not 
known. In addition to water shedding, ice shedding is a 
major consideration for situations involving high speeds, 
rotorcraft and deicing systems. Olsen and Newton 8 and 
Ruff 6 have discussed the similitude requirements for 
shedding. 

Both of these shedding phenomena take place away 
from stagnation, which is where the scaling analyses are 
focused. Thus, the present approach to scaling cannot 
easily incorporate them. However, for some conditions 
these events might need to be considered for accurate 
simulation of horn shape. Until high-speed, close-up 
studies can better identify the processes taking place, 
shedding can not be analyzed. Fortunately, for typical 
Appendix C conditions, reasonable success with scaling 
has been demonstrated without including a shedding 
analysis. 

CONCLUSION 

This paper reviewed the similarity parameters identified 
by various research studies from the standpoint of icing 
scaling. Scaling methods incorporate the physics of ic- 
ing by requiring that scale and reference values of the 
most important similarity parameters should match. 

From the drop momentum equation, Langmuir and 
Blodgett 7 derived the modified inertia parameter, K 0 , 
given by equation (3). The stagnation collection effi- 
ciency at the leading edge, p 0 , is directly related to K 0 by 
equation (8). Main ice shapes are not particularly sensi- 
tive to Po, so some flexibility in matching this parameter 
appears to be possible. By matching scale and refer- 
ence values of K 0 or p 0 the scale drop size can be de- 
termined. 

From an analysis of the water catch the accumulation 
parameter can be formulated, and from it the scale icing 
time can be determined (equation (12)). 
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Messinger’s 12 freezing fraction, n 0 , of equation (28) was 
derived from the surface energy balance. n 0 has a 
strong effect on the glaze ice horn angle and the lead- 
ing-edge thickness. A second parameter from the en- 
ergy balance is </>, the water-energy transfer parameter, 
given by equation (26). The Ruff scaling method 5 6 7 
matches <f> s to </) R to find scale temperature and n 0 ,s to 
n 0 , R to find scale LWC. However, ice shapes appear to 
be much less sensitive to <p than to n 0 . The expression 
for freezing fraction presented in this paper was vali- 
dated by comparing analytical and experimental values. 
The appearance of experimental ice shapes was also 
consistent with the analytical values of freezing fraction: 
when the freezing fraction was unity, ice shapes had the 
typical took of rime ice, and when n 0 was less than unity, 
glaze features appeared at the leading edge. 

Consideration of surface-water dynamics led to the con- 
clusion that an additional similarity parameter is needed 
to find scale velocity. A Weber number based on free- 
stream velocity and model size (eq. (32)) was presented, 
but there are other possible Weber numbers as well as 
parameters such as the non-dimensional water-film 
thickness that could be used. Because there is not a 
significant difference between the values of velocity 
found using most of the proposed parameters, it is not 
possible to discriminate between them by evaluating 
scaling results. Additional studies of surface phenom- 
ena, including water-film behavior and drop impact and 
splashing, are needed to understand the physics better. 
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NOMENCLATURE 


A c accumulation parameter, dimensionless 

AOA angle of attack, ° 

b relative heat factor, dimensionless 

c airfoil chord, in 

c p constant-pressure specific heat, Btu/lbm R 

C D drag coefficient of drop, dimensionless 

d cylinder diameter or twice the airfoil leading- 
edge radius, in 

D v diffusivity of water vapor in air, ft 2 /s 
e 1, e 2, e 3 powers in equation (30), dimensionless 
h c convective heat transfer coefficient, Btu/hr ft 2 R 
hfiim water-film thickness, in 

h G gas-phase mass transfer coefficient, Ibm/hr ft 2 

k thermal conductivity, Btu/hr ft R 

K drop inertia parameter, dimensionless 

K 0 modified inertia parameter, dimensionless 

L undefined length proportional to model diameter 

or chord, in 

LWC liquid water content of cloud, lbm/ft 3 

m mass flux of water per unit time, lbm/ft 2 s 

rh e mass flux of water evaporated, lbm/ft 2 s 

M Mach number, dimensionless 

n freezing fraction, dimensionless 

n a analytical value of freezing fraction, dimen- 

sionless 

n e freezing fraction found from leading-edge thick- 

ness of experimental ice shapes, dimensionless 
Nu Nusselt number, dimensionless 

p pressure, psi 

p ww vapor pressure of water over liquid water, psia 
p„ vapor pressure of water in the atmosphere, psia 
Pr Prandtl number of air, dimensionless 

q c surface heat loss due to convection, Btu/hr ft 2 
q e surface heat loss from evaporation, Btu/hr ft 2 
q f surface heat gain from release of latent heat of 
fusion, Btu/hr ft 2 

q k surface heat gain from kinetic energy of water 
drops, Btu/hr ft 2 

q w surface heat loss to raise temperature of imping- 
ing water to freezing, Btu/hr ft 2 
r recovery factor, dimensionless 

r, e leading-edge radius, in 

r A evaporative-to-convective heat transfer term, 
dimensionless 

R a gas constant for air, Ibf ft/lbm R 

Re Reynolds number, dimensionless 

Re re i Reynolds number of drop relative to airstream, 
dimensionless 


Re s Reynolds number of drop, dimensionless 
s Surface distance from airfoil stagnation, in 

Sc Schmidt number of air, dimensionless 

t temperature, °F 

t b/ temperature of air in boundary layer, °F 

t f freezing temperature of water, 0°F 

tfiim film temperature, °F 

T absolute temperature, R 

T f absolute freezing temperature of water, R 

u air velocity vector, dimensionless 

V free-stream velocity (air speed), mph 

\Ne Weber number based on drop diameter, dimen- 
sionless 

We c Weber number based on model length dimen- 
sion, dimensionless 

We L Weber number based on length L, dimen- 
sionless 

x horizontal coordinate, in 

y vertical coordinate, in 

P catch efficiency, dimensionless 

y ratio of specific heats for air, 1 .4 

5 drop median volume diameter (MVD), pm 
A ice thickness, in 

0 air energy transfer parameter, R 

A drop range in absence of gravity, ft 
^States drop range in absence of gravity if Stokes' law 
applies, ft 

A f latent heat of freezing, Btu/lbm 

A v latent heat of vaporization, Btu/lbm 
p viscosity, lbm/ft s 

p density, lbm/ft 3 

u wa Surface tension of water over air, Ibf/ft 

t icing time, min 

<j) drop energy transfer parameter, R 

General Subscripts: 

0 stagnation value 

a air 

/ ice 

R reference conditions 

s surface 

S scale conditions 

st static 

tot total 

w water 
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